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Abstract. The ALICE experiment has measured the D meson production in pp 
and Pb-Pb collisions at the LHC at y/s = 7 and 2.76 TeV and y/sj^ =2.76 TeV 
respectively, via the exclusive reconstruction of hadronic decay channels. The analyses 
of the D° — > K~7r + and D + — > K~tt + tt + channels will be described and the preliminary 
results for the D° and D + nuclear modification factor will be presented. 



1. Introduction 

The comparison of heavy flavour production in proton-proton and heavy-ion collisions 
allows to probe the properties of the high-density QCD medium formed in the latter 
and to study the mechanism of in- medium partonic energy loss [TJ. A sensitive 
observable is the nuclear modification factor, defined, for a particle species h, as 
^AA(Pt) = (TAA^xdi-^/dpt ' wnere -^aa * s the normalized yield measured in heavy- ion 
collisions, (Xaa) is the average nuclear overlap function calculated with the Glauber 
model in the considered centrality range, and a^ p is the production cross-section in 
pp collisions. For hard processes, in the absence of any medium and cold nuclear 
matter effect, -Raa = 1 is expected. By comparing the nuclear modification factors 
of charged pions (-Raa)' mostly originating from gluon fragmentation, with that of 
hadrons with charm (-Raa) an d beauty (-Raa) the dependence of the energy loss on 
the parton nature (quark/gluon) and mass can be investigated. A mass ordering 
pattern i? AA (p t ) < i?° A (p t ) < -RaaG°0 has been predicted [2]. In these proceedings the 
measurement of the D° and D + -Raa in Pb-Pb collisions at the LHC is presented. More 
information on the ALICE detector can be found in [3j H] ; the procedure to determine 
the collision centrality is described in [5]. 

2. Measurement of D meson production with the ALICE detector 

The production of D° and D + (cr ^ 123 and 312 /mi respectively [6]) was measured in 
pp and Pb-Pb collisions at central rapidity (\y\ < 0.5) via the exclusive reconstruction of 
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the decays D° -> K-vr+ (with branching ratio, BR=3.89±0.05% 0) and D+ -> K-7T+7T+ 
(BR=9.4 ± 0.4% [6]). The analysis strategy for the extraction of the D° and D + signals 
out of the large combinatorial background from uncorrelated tracks is based on the 
reconstruction and selection of secondary vertex topologies with significant separation 
(typically a few hundred micrometer) from the primary vertex. The Time Projection 
Chamber (TPC) and the Inner Tracking System (ITS) detectors provide a spatial 
resolution on the track position in the vicinity of the primary vertex of the order of few 
tens of microns at sufficiently high p t [7] . Tracks displaced from the primary vertex are 
selected to reconstruct D meson candidates. Good alignment between the reconstructed 
meson momentum and the flight direction between the primary and secondary vertex 
is also required. The identification of the charged kaon in the TPC and Time Of 
Flight detectors provides additional background rejection. A particle identification 
strategy that preserves most of the D meson signal was adopted. Similar analyses 
were performed on pp and Pb-Pb data, with a tighter selection in the latter case, 
dictated by the higher combinatorial background. To extract the signal, a fit to the 
invariant mass distribution is performed, using a Gaussian function for the signal peak 
and an exponential function for the background. Figure Q] (left) shows the invariant 
mass distribution of D° candidates for p t > 2 GeV/c, obtained from the analysis of 
ps 3 x 10 6 Pb-Pb events in the — 20% centrality range. The 'raw' signal is corrected 
for acceptance and efficiency using Monte Carlo simulations based on Pythia (Perugia-0 
tuning) and HIJING event generators. The contribution of D mesons from B decays was 
evaluated relying on the FONLL prediction, which describes well bottom production at 
the Tevatron [8] and at the LHC [10]. In the Pb-Pb case, the FONLL prediction of 
secondary D mesons production in pp collisions is multiplied by (N co \\) x R® A , where 
the hypothesis on the B meson Raa encodes all potential nuclear and medium effects 
affecting B production. The systematic error due to the FONLL theoretical uncertainty 
on beauty prediction was estimated from the spread of the results recovered using the 
minimum and maximum predictions for secondary D production. As reported in [11] , the 
D meson production cross-sections in pp collisions at yfs = 7 TeV are well described by 
predictions based on pQCD calculations as FONLL [8] and GM-VFNS [12]. To provide 
the reference cross- section at 2.76 TeV, needed to compute the -Raa, the measurement at 
7 TeV is scaled by the ratio of FONLL predictions for D meson production at \fs = 2.76 
and 7 TeV [13J. The uncertainty on the scaling decreases with p t from 25% to 10%, 
as estimated by varying the parameters entering the FONLL calculations, namely the 
factorization and normalization scales and the charm quark mass. In March 2011 a 
sample of ~ 6.5 x 10 7 events from pp collisions at y/s = 2.76 TeV was collected. The D° 
and D + signals were measured in the 2 < p t < 8 GeV/c range. While the accumulated 
statistics did not allow for determining a pp reference over the whole momentum range, it 
served as an important cross check of the theoretical scaling procedure in the momentum 
range where the data sets overlap. Considering only the statistical errors, the ratio 
between the cross-section measured at 2.76 TeV and that obtained from the scaling of 
the measurement at 7 TeV is compatible with 1 within la for the D + while, for the 
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Figure 1. Left: invariant mass distributions of D° candidates for pt > 2 GeV/c 
obtained from the analysis of ~ 3 x 10 6 central (0 — 20%) Pb-Pb events at v^nn = 
2.76 TeV. The curves show the fit functions as described in the text. Right: D° spectra 



for central (0 — 20%) and peripheral (40 



Pb-Pb events at Vs NN = 2.76 TeV 



compared to the reference spectra obtained by scaling the ALICE measurement in pp 
collisions at 7 TeV (see text). 



D°, within la above 4 GeV/c and within 2a for 2 < p t < 4 GeV/c. A further 20% 
asymmetric contribution to the systematic error in the 2 < p t < 4 GeV/c range has 
been added on the D° cross-section at 2.76 TeV used as the pp reference in the Raa 
calculations. Figure [1] (right) shows the D° p t spectra measured in central (0 — 20%) and 
peripheral (40 — 80%) Pb-Pb collisions, compared to their respective reference spectra. 

3. D meson nuclear modification factors 

The D meson production is suppressed by a factor 4-5 in central events for p t > 5 GeV/c, 
as quantified by the nuclear modification factors shown in Fig. [2] (left). The D° and 
D + Raa agree within errors. For the D°, the statistical uncertainty is of the order of 
20 — 25%, the estimated total systematic uncertainty, accounting for the uncertainties on 
the signal extraction procedure, on track reconstruction efficiency, on the MC corrections 
for reconstruction acceptance, cut and PID selection, and for possible differences in 
the D° and D° reconstruction varies from ^26% for 2 < p t < 4 GeV/c to ±25% 
for 8 < p t < 12 GeV/c. The hypothesis on i?^ A is varied in order to span the 
range 0.3 < -Raa/-^aa < 3, with the -Raa/^aa calculated a posteriori. The range 
of R^aa vaiues obtained in each p t bin is considered as the systematic error due to 
the -Raa assumption. For the D°, the maximum spread is of the order of 15% in 
the bin 2 < p t < 4 GeV/c, 25% in the bin 4 < p t < 5 GeV/c and 10% in the bin 
8 < p t < 12 GeV/c. 

In central collisions, -Raa decreases with p t from ~ 0.45 at low p t to w 0.2 at high 
p t . The D meson -Raa is compatible within errors with the charged pion -Raa: the larger 



4 



1.8- 
1.6- 



Pb-Pb Vs^=2.76TeV 



ALICE Preliminary 



* D u R AA 0-20% CC " 
■ D + R AA 0-20% CC " 

• rc* R.. 0-20% CC - 




1.2 



12 14 

p [GeV/c] 



Q. 

E 
o 



?0.8 
0.6 
0.4 
0.2 



ALICE Preliminary 



■ R AA and statistical unc. 
| Data systematic unc. 
j H Correlated systematic unc. 



" Pb-Pb,\/s N11 =2.76 TeV, 6 < p T D° < 12 GeV/c 



60-80% 



40-60% 



20-40% 



10-20% 0-10% 

Centrality 



Figure 2. Left: D° and D + nuclear modification factor measured in central (0 — 20%) 
Pb-Pb collisions. The vertical lines represent the statistical errors. The filled grey 
box around Raa = 1 represents the uncertainty on the normalization of the pp cross- 
section and on (Taa)- The empty boxes represent the total systematic uncertainties. 
The charged pions Raa is shown for comparison. Right: D° nuclear modification factor 
as a function of centrality for 6 < < 12 GeV/c. 



statistics expected from the higher luminosity Pb-Pb 2011 run and a reduction of the 
systematic errors should allow to conclude whether R® A > -Raa f° r Pt ~ 5 GeV/c 
and quantify the difference. The dependence on the centrality is shown in Fig. [2] 
(right) for 6 < p t < 12 GeV/c: Raa decreases from rs 0.7 in peripheral (60 — 80%) 
to ~ 0.2 in central (0 — 10%) events. As verified by calculating the Raa expected 
from pQCD calculation based on the MNR code [H] and nuclear PDF from EPS09 
parametrization [15], nuclear shadowing yields a relatively small effect for p t > 5 GeV/c 
(0.85 < i?° A < 1.08 at p t = 5 GeV/c). Therefore, the observed suppression can be 
considered as an evidence of in-medium charm quark energy loss. 
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